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ABSTRACT 

We calculate the dependence of the electrostatic persistence length, l e , of weakly charged 
flexible polyelectrolyte chains using a self-consistent variational theory. The variation of l e 
with k, the inverse Debye screening length, is controlled by the parameter IqIb/A 2 where Iq 
is the bare persistence length, Ib is the Bjerrum length, and A is the mean distance between 
charges along the chain. Several distinct regimes for the dependence of l e on k emerge 
depending on the value of IqIb/A 2 . We show that when IqIb/ A 2 <C 1 we recover the classical 
result, l e oc k~ 2 . For intermediate values of IqIb/A 2 , l e oc In this regime one can also 
get l e oc k~ v with y < 1 depending on the strength of the Coulomb interaction. Qualitative 
comparisons between our theory and simulations as well as other theories are presented. 
PACS number(s): 36.20.-r, 33.15.Bh, 64.60.-i, 41.20.-q 
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1 Introduction 



In recent years, significant attention has been paid to the old problem of the ionic dependence 
of the size of isolated polyelectrolyte chains [1-4]. The size of the chain is most conveniently 
expressed in terms of the electrostatic persistence length. The computation of the persistence 
length of a polyelectrolyte chain is difficult because of the interplay of several length scales. 
Significant steps in the computation of the size of charged polymers was taken by Odjik || 
and Skolnick and Fixman || (OSF), who introduced the concept of the electrostatic persis- 
tence length. These authors considered a stiff polyelectrolyte chain near the rod limit so that 
only small excursions from linear conformations are relevant. By assuming that charges on 
the rod- like polyelectrolyte interact via the Debye Hiickel potential, = /se _Kr /r, where 
Ib = e 2 / A-KeoksT is the Bjerrum length, A is the average distance between charges, k _1 is the 
screening length, OSF showed that the total persistence length of the chains can be written 
as l p = lo + l e - Here l is the bare persistence length, and the electrostatic persistence length 
l e in the OSF theory is 

l OSF = \^ (1.1) 

It is believed that as long as the chain is stiff the dependence of l e on k given above is correct 
[2] . The above result is exact in large dimensions regardless of the intrinsic stiffness of the 
chain backbone. 

Many experiments on polyelectrolytes have been interpreted using the OSF theory. How- 
ever, light scattering experiments on some polyelectrolytes systems do show deviations from 
the OSF theory |7|, |8|, ||. The major approximation in the OSF treatment is that the poly- 
electrolyte is rod-like, and hence deviations from linear configuration is small. If the intrinsic 
backbone of the polymer chain is flexible, so that large excursions from linearity are possi- 
ble, it is likely that the OSF result may not be valid. In fact, the validity of Eq. (1.1) for 
chains that are intrinsically flexible has been questioned using theoretical arguments [7-9] 



and computer simulations [13|. 14]. Starting with the work of Barrat and Joanny [[l l| (BJ), 



several variational calculations fllO| , |l"2"f show that, when the chains are intrinsically flexible, 
the electrostatic persistence length l e ~ k -1 . We had shown earlier |10j that the predicted 
is valid when the controlling parameter IqIb/A 2 ~ 0(1), i.e., when the chain is in the 



so called non-asymptotic regime |10fl . Li and Witten [TJ[ argued that, unlike the variational 



theories |IJ |TT|, |12| ) approximate inclusion of fluctuations still leads to the result predicted 
by OSF even if the backbone is flexible. In other words, the OSF result given in Eq. (1.1) 
is always valid. In the above referenced variational theories it is generally believed that the 
electrostatic interactions generally stiffen the chain significantly so that the trial Hamilto- 
nian should consist of terms that account for the electrostatic-mediated interaction rigidity. 
In contrast, Muthukumar has examined the influence of electrostatic interactions on the 



size of the flexible polyelectrolyte chains [[U| • This theory, which also emphasize the role of 
excluded volume interactions, has been used to interpret experimental data [[Tj|. A precise 
form of the ionic dependence of the electrostatic persistence length that covers all values of 
stiffness has not emerged from any these studies. 

More recently, Micka and Kremer [13] have extracted the persistence length of weakly 
charged polyelectrolytes of varying stiffness using Monte Carlo simulations. They argue that 
there is no apparent simple scaling of l e with k in all regimes. This raises the possibility that 
a simple concept of electrostatic persistence length may not be useful in determining the size 
of polyelectrolyte chains under diverse conditions. The numerical results can be apparently 
be rationalized using a field-theoretical renormalization group analysis of flexible chains. 
Liverpool and Stapper JL7] suggest that the scaling form for the mean square end-to-end 
distance, (-R 2 ), provides a good fit of the computer simulation data. 

The preceding survey suggests that we still do not have a simple theory for computing l e 
for charged polymer chains with arbitrary stiffness. The purpose of this paper is to reexamine 
the variational type theories for computing l e using the trial Hamiltonian suggested by BJ 
for flexible chains |TTJ. Rather than compute the approximate free energy using the Gibbs- 
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Bogoliubov inequality as BJ did, we calculate (R 2 ) directly using the uniform expansion 
method HT9| . We show explicitly that when l l B /A 2 <C 1, the OSF result given in Eq. 
(1.1) is recovered. The rest of the paper is organized as follows. In Sec. 2, we present 
the variational theory for the electrostatic persistence length in the asymptotic limit, i.e., 
when IoIb/A 2 <C 1. The case when IqIb ~ C(i)> which is most appropriate for examining 
the crossover regime, is presented in Sec. 3. Comparisons of the theoretical predictions to 
the simulation result are discussed in Sec. 4. The paper is concluded in Sec. 5 with a few 
additional remarks. 



2 Variational Theory: Chain under tension 



The flexible polyelectrolyte chain consists of iV segments each of length Iq on an average. 
The charges on the chain are assumed to be a distance A apart. The interaction between the 
charges on the chain is taken to be given by the Debye-Hiickel (DH) potential. In this model 
the counterions are not explicitly treated. In the continuum limit the partition function of 
the polyelectrolyte chain is taken to be 



V[t(s)]W[t(s) 



(2.1a) 



where the statistical weight associated with a given realization of the polyion described in 
terms of the curve r(s) is 



W[r(s)] = W°[r(s)]exp 



44 /Yd,d^ 



-re|r(s')— r(s)| 



2 A 2 Jo Jo r(.s') - r(.s) 

and W°[r(s)] is the weight for the chain backbone in the absence of charges 



(2.16) 



W°[r(s)] =exp 



3 r L ( dr(s) 



(2 Ac) 



21 Jo V ds J 

with I = 2l . In order to eliminate the self-interaction in the DH potential we impose the 
condition that \s' — s\ > so where so ~ Iq. The cut-off dependence vanishes in the variational 
theory described below. 



The direct evaluation of Z in Eq. (2.1a) is not possible due to the non-Markovian nature 
of the electrostatic interaction term. Consequently, we have to resort to an approximate 
calculation. Here we use a variational method, first proposed by Edwards and Singh fl^| 
(ES), in their treatment of the standard excluded volume problem in neutral flexible poly- 
mers. The basic idea of the ES method is to replace the weight by a simpler trial weight 
with undetermined parameters. The parameters are chosen so that the size (say the mean 
square end-to-end distance) calculated using the trial weight, (-R 2 ) t becomes identical to 
that given by the true weight, (-R 2 ). The practical implementation of this method requires 
that the interaction term can somehow be treated as small so that the deviation between 
(R 2 ) t and (R 2 ) can be computed in powers of the small parameter. The ES method has 
been successfully used to predict the conformational properties of flexible polyions under a 



variety of solvent conditions [16, 21 



A word of caution concerning the ES method (and in general variational methods) is in 
order. Firstly, the results of variational calculations are only as good as the chosen trial 
Hamiltonian. In polymer problems there is an additional distinction. ES advocate that the 
variational parameters be chosen to ensure (-R 2 ) t = (-R 2 ) rather than by employing the Gibbs- 
Bogoliubov inequality, i.e., by minimizing the free energy. For polymers with short range 
interactions it has been shown that the ES method is superior to others that minimize the 
free energy by a suitable variational calculation. For problems with long range interactions 
the direct minimization of free energy can give reliable results for the scaling of (R 2 ) with 
N. For the present problem the range of interaction is neither short nor too large. From 
this discussion it should be clear that the utility of the variational methods should be viewed 
with some caution. 

The choice of the trial weight or Hamiltonian is extremely important. An appropriate 
trial Hamiltonian is chosen by considering the case of k — 0. Here, the weakly charged 
flexible polyelectrolyte can be considered as an array of blobs of size D. The structure of the 
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chain inside a blob resembles that of a non- interacting (i.e., the electrostatic interactions are 
unimportant) while there is ordering among the blobs due to the Coulomb repulsion beyond 
the length scale D. The physical picture of the charged chain is that of one under tension Q, 
where the blob size (short length scale effect) is determined by the magnitude of the (as yet 
undetermined value of the tension), while at long length scales the chain is a free chain 
with an effective step length different from Iq. Thus, in order to determine the structure 
of such polyions, it is crucial to incorporate the presence of two length scales, namely, D 
that describes the local ordering and correlations between the semiflexible blobs that gives 
the effective persistence length. Based on this physically appealing picture, Barrat and 
Joanny |TT| proposed a trial probability weight W t [r(s)] associated with the configuration 
[r(s)] that succinctly captures the conformation of weakly charged polyions described above. 
In particular they chose 



W t [r(s)] = W°[r(s)] / D[t(s)]W[t(s)]exp 



rl L dst(s) 9r(s) 



o 



(2.2) 



ds 

where W°[r(s)] is the Boltzmann weight for the non-interacting chain (Cf Eq. (2.1c)), t(s) 
is the unit vector that gives the direction of the tension at position s. The weight Wt[r(s)] 
is associated with a particular realization of the direction of the chain tension and is given 
by 



W[t(*)]= J] 5(t( S ) 2 - 1) exp 

0<s<L 



(2.3) 



2 Jo V ds J 

The weight associated with the tension vector t(s) is an attempt to capture the physics that 
the charged chain on a scale rj behaves as one under tension. On the scale rj the chain is 
adequately stiff so that it can be described as worm-like for which the Boltzmann weight 
is approximately given by Eq. (2.3). The finite range of correlations among the blobs is 
described by the correlation in t(s) such that 

(t(s)-t(s , )) w = exp[-\s'-s\/r ] ]. (2.4) 

The trial weight has two parameters namely rj (the effective persistence length) and the chain 



tension r. Both these are adjusted so that the deviation between (R 2 ) and {R 2 ) t is zero to 
first order in the interaction (second term in Eq. (2.1b)) term. 

From the form of the trial weight W t [r(s)] it is clear that the blob size D ~ r _1 . This 
can be formally established by computing (|r(s') — r(s)| 2 ) Wt which can be obtained from the 
generating functional defined by 

£(k;s',s) = ( e ik -( r ( s ')-rW)\ _ (2 .5) 

\ / Wt 

The evaluation of £(k; s', s) up to order k 2 is given by 

£(k;s',s) = e -^'^ / e ^/3-/;'t( S ")d,"\ 

\ / w t 

= _ i k H\ s > _ s | _ fj* d Sl ds 2 (k • t( Sl ) k • t(s 2 )) m 

= _i A ; 2 /| S / - S |-^A; 2 / 2 r 2 ff d Sl ds 2 (t( Sl ) • t(s 2 )) m . (2.6) 

J S J S 

The mean square internal distance is 

(\r(s')-r(s)\ 2 ) m = -Q(d/dk 2 )\ k=0 C(k;s',s) 

= \ s > _ s | + iZ 2 r 2 [277|s' - s| - 2^ 2 (1 - e -l a '- a l/")], (2.7) 

where we have used Eq. (2.4). By equating the two terms in Eq. (2.7) we can infer the 
contour length S ~ 1/Zt 2 and the blob size D ~ (SI) 1 / 2 ~ r _1 . The blob size is naturally 
expressed also in terms of Iq, /b, and A and is given by Eq. (2.7). For long chains i.e., L ^> S 
we can neglect the first term in Eq. (3.7), and the size of the chain {R 2 ) Wt becomes 

(R 2 ) Wt = ML) - r(0)| 2 ) Wt = 2l p L' - 2l 2 p (l - e"^), ( 2 .8) 

where l p ~ r/lr is the physically observable persistence length and L' is the longitudinal size 
of the corresponding stretched chain 

L' w W~ X L w (WbM 2 ) 1/3 I. (2.9) 
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Note that the notion of stretched blobs is only valid for l v /D >> 1. 

In order to set up the calculation of l v using the ES method we write 

(R 2 ) = (R 2 ) m + AR 2 (l p , Z B , A, k) (2.10) 

where (R 2 ) is to be computed using the weight given in Eq. (2.1). As explained before the 
variational parameters 77 and r are determined from the condition AR 2 = 0. To obtain an 
explicit expression for AR 2 let us write 



JgjrCg)] it> 2 W[r(s)] 
^ ; jT>[r( S )] W[r( S )] 

JX>[r(s)] W t [r(g)] # 2 exp{-£[r(s)]} 
/P[r(s)] W t [r(s)] exp{-B[r(s)]} 

where 



(2.11a) 



/„ /■£ f L e -«|r(s')-r(s)| ,L „ ,L //}t\ 2 

If W t [r(s)] is a good approximation to W[r(s)] then B[r(s)] can be treated as small. Conse- 
quently up to 0(B[r(s)]) the vanishing of AR 2 leads to the condition: 

(R 2 B[r(s)]) Wt = (R 2 ) m (B[r(s)]) Wt . (2.12) 

By using the following equations 



w t \ ^° / Wt 



eq. (2.12) may be written as 

'p J-Wv* " ^<^>w t = (^>w t - (^W^V (2-14) 



where V is the total interaction energy given by the first term in Eq. (2.11b). Using Eq. 
(2.8) for {R 2 ) m the L.H.S. of Eq. (2.14) can be easily calculated. 



The computation of {R 2 V) Wt and (V) w is somewhat involved. Here we provide of a 
sketch of the calculation. The details are given in Appendix A. An expression of (R 2 V) Wt 
can be obtained in terms of the generating functional 



£i(k;0,L|q; s, s') = ( e *W)-*(.o))+i*W)-'W) 



(2.15) 



with the result 

(R 2 V) 
Similarly 



u 



dsds' 



-6- 



d 



d 3 q Air 

(2tt) 3 q 2 + k 2 L "<9k 2 



k=0 



£(k; 0, L|q; s, s'). 



(2.16) 



-6 



g 
dk 2 



£(k; o, L)^ // d sds - / ^j-^Aq; «') (2.17) 



- (aov + K 2 

where £(q; s, s') is given by Eq. (2.5). 

In order to obtain a tractable expression for l p we make an additional, physically mo- 
tivated, approximation to evaluate Eqs. (2.16) and (2.17). On a scale larger than rj the 
chain is effectively Gaussian provided L/rj 3> 1. Thus, the effect of stiffness, arising from 
stretching due to electrostatic repulsion, is important only on the scale of the order rj. The 
chain on the scale of rj can be treated as essentially rod-like with minimal conformational 
fluctuations. In order to obtain the generating functional C and C\ we need to compute 
terms with the structure (V Q '^o • With the approximation that on scale of the order 

of \s' — s\ = i] the chain behaves essentially as rod-like we can write 



iQ-J^t(s)ds 



w t 



e ^ Js 

Wt \ / rod 



(2.18) 



With this physically motivated approximation the generating functional can be calculated. 
The result for the persistence length can be written as (after considerable algebra) 

d 3 g 



* %h II dsds 'I ^ e ~ q2ll " sl/e {^ s> - *) V q " /3 - /;V)d 0' 



d 



k=0 



(k+q)lT/3.// t(s")ds"\' _ /-iklT/3-/; t( S ")d S "\V iqZr/3-/; t(s")ds"\' 



(2.19) 



where we have neglected numerical factors. The prime in Eq. (2.19) indicates that the 
statistical average (V Q '^o t<s ^ ds ) w has been performed using the factorization described in 
Eq. (2.18). The evaluations of various quantities in Eq. (2.19) are given in Appendix A. 
We now provide an argument that only the quantity arising from the local interaction, i.e., 



length. The long length scale term provides an interaction of the excluded volume type and 
is not relevant in the calculation of l p . In order to show this let us split the integrals in 
Eq. (3.19) into those arising from \s' — s\ < rj and into the contributions arising from scales 
larger than rj. The integrals involving the tension term in t(s) are modulated by the Gaussian 
factor exp(— q 2 l\s' — sj/6) which accounts for local fluctuations inside a blob. First consider 
the limit qlr\s' — s\ <^1. In this limit we can approximate 



for \s' — s\ < rj. It is clear from Eq. (2.20) that on scales less than S (~ I/It 2 ) the Gaussian 
term is more important (see also Eq. (2.7)) while on the scale S < \s' — s\ < i] Eq. (2.20) 
dominates. On this length scale the segment of the chain behaves as a rod, and hence charges 
separated by distance \s' — s\ = 0(rf) can be assumed to interact via the screened Coulomb 
interaction. Now consider scales larger than rj. On these long length scales the chain can 
be pictured as a sequence of rods of length l p ~ r/lr surrounded by an impenetrable charge 
cloud of diameter d. Two such segments can come arbitrarily close together to a distance of 
the order of d because the chain is Gaussian on long length scale. The interaction of such rod 
segments gives rise to a short range excluded volume type terms. Since these interactions are 
not relevant for describing the local ordering (on scale 77) we only include contributions from 
short length scales in Eq. (2.19). Polyions on length scale greater than r/ can be considered 
as flexible chains with an effective step length l p interacting via an effective excluded volume 
type interactions. 

If we substitute Eqs. (A. 5), (A.7), (A.8), and (A.ll) into Eq. (2.29) we obtain the result 





J / Wt 




(2.20) 
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for l p , which is valid for the screened polyions in the asymptotic limit defined as IqIb <C A 2 , 

<>~(5D 2 |( fl +^). < 2 - 21 > 

where (S/A) 2 (l B /D) ~ 0(1) is the electrostatic energy of a blob. The second term in Eq. 
(2.21), which varies as k~ 2 , can be identified as the electrostatic persistence length while the 
^-independent term 0(1)D can be identified as a renormalized bare persistence length due 
to the chain fluctuations. Surprisingly, this result is quite similar to the one derived from 
the OSF result for l p by replacing l , A by the blob size D and the charge per monomer e by 
the charge ge of a blob, as was done by Khokhlov and Khachaturian [^] (KK) without the 
benefit of any explicit derivation: 

lp - lo + 4^A 2 D+ {a) D D^ 2 

= £> + 0(l)_L (2.22) 

A derivation of the result displayed in Eq. (2.22) was also recently provided by Li and Witten 
[15] using an entirely different approach. From the preceding calculations, we find that the 
fluctuations in chain conformations of asymptotic polyions (IqIb/A 2 <C 1) do not invalidate 
the qualitative aspect of the OSF result. They effectively reduce the direct distances between 
charges so as to renormalize the bare parameters Iq, A into D and e into ge, as was first 
recognized by KK. In this limit, local fluctuations inside the blobs are strong enough to reduce 
significantly the direct distance between two consecutive charges by the factor (1 Ib/ A 2 ) 1 / 3 
which is much less than 1. This results in much stronger effective Coulomb repulsion between 
two consecutive charges. Since the length scale of l p is much larger than the blob size D and 
ordering as implied in these cases does not refer to the local structure of blobs, we can expect 
that the local fluctuations inside a blob not to affect significantly the property of the chain 
much beyond it. As a result we obtain qualitatively similar dependence of l p on k as that 
obtained for stiff polyions. Note that this is relevant only for the asymptotic case where each 
blob contains a large number of segments which roughly obey Gaussian statistics. Indeed 
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computer simulations performed by Higgs and Orland [E^] correspond to the asymptotic 
case (IqIb/A 2 = 0.05), showing that l p is much larger than k" 1 . Their results support the 
expression for l p given in Eq. (2.22). 

The persistence length of unscreened polyions, however, diverges as L — > oo. This is 
because all charges on the chain contribute to the chain stiffening. In the OSF theory it 
is proportional to 1-qL 2 /A 2 . To obtain the corresponding persistence length for unscreened 
flexible polyions, we replace Ib/A 2 by 0(1)/D and L by V 

r/2 

lp ~ -p- (2-23) 

It is interesting to contrast the results given in Eq. (2.22) with the calculations of BJ. Our 
theory and the BJ analysis use exactly the same variational Hamiltonian. The difference is 
that BJ identified the persistence length by obtaining an upper bound to the free energy of 
the system using the Gibbs-Bogoliubov inequality. We, on the other hand, have computed 
(R 2 ) using the ES method. The persistence length of the chain is a "local" measure of the 
size of the polymer. In the formulation of BJ, contributions from these scales are averaged 
out whereas the ES method, at least approximately, accounts for short length scales effects. 
We believe this is the major source for differing the conclusions drawn by BJ. This point is 
further illustrated in Appendix B where we show the deficiencies of the variational theory, 
using the free energy for computing "local" quantities such as persistence length. 



3 Polyions with WbM 2 ~ 0(1); Gaussian Chain Model 

In the previous section we considered the case of flexible polyelectrolyte in the so called 
asymptotic limit i.e., when A 2 3> Zo^b- If the magnitude of A, l , 1-q are similar, then one is in 
the so called non-asymptotic limit and represents crossover regime from small values to large 
values of IqI-q/A 2 . In the limit of IqI^/A 2 ^> 1, which is the best defined limit, the electrostatic 
persistence length is given by the Odjik-Skolnick-Fixman theory. When IqL-q/A 2 ~ 1 then it 
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is clear that (in all likelihood) the blob size, D, is of the same order as Iq or D/l ~ 1. In this 
limit the distinction between the blob size and the bare persistence length is meaningless 
and the chain becomes ordered on scales greater than l due to electrostatic repulsion. As 
a result it seems reasonable that one can choose a simpler trial weight in contrast to the 
polyions in the asymptotic limit which treated in the previous section. 

Following the arguments given above, we chose the following trial weight for IqI-q/ A? ~ 1 



Wi = exp 



3 

2h 



'dr(s) 
ds 



(3.1) 



which is a Gaussian (non-interacting) theory with an effective step length / x . The self- 
consistent equation using the ES procedure, namely one that gives vanishing AR 2 to first 
order in 



7 - - - K |r( S ')-r(s)| 

B[r(s)] = ^ I I dsds' 



2A 2 



\r(s') — r(s) 



K)/:©-< <"> 



for 1 1 is 



{R 2 V) m - {tf) m {V) 



(3.3) 



The calculation of the various quantities using the Gaussian weight is straightforward, and 
using these results Eq. (3.3) becomes 

1 1 1 Zb /■/",,,/ , „ f d 3 g Anq 2 



II 



dsds'is' 



-q 2 h\s'~s\/& 



(3.4) 



/ h 18LA 2 JJ y ' J (2vr) 3 q 2 + k 2 

The evaluation of the integral in Eq. (3.4) is done by imposing a lower cut-off for q = q = /j" 1 . 
For this choice of go the contribution from charges separated by distance larger than l\ is 
exponentially small. Furthermore, since li is the only scale in the Gaussian weight it is 
reasonable to set q = Zf 1 . With this cut-off the evaluation of the integrals in Eq. (3.4) is 
straightforward and one obtains the following equation for l\ 



I 2 



1 1 



ah 



.1 IJ A 2 K 2 
where a is a numerical constant. 



1 



IT 



tan 



1 

1 X K 



(3.5) 
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The self-consistent equation can be easily solved if l\ is somewhat larger than the Debye 

screening length n~ l . In this case, we can neglect the underlined term in Eq. (3.5). This 

results in the following solution for l p = l\j1 

lo + ^osf, if ^osf -C lo 

*P~< . . (3-6) 



(l l B /A 2 ) 1/2 K -\ if Iqsf^I 







where Zosf ~ Ib/A 2 k 2 is the electrostatic persistence length of stiff chains derived from 
the variational equation. For the flexible chains (/osf ^ or equivalently li ^> l ), the 
persistence length varies as k , which is in accord with Barrat and Joanny [13] and our 
earlier calculation [10]. The result in Eq. (3.6) has also been obtained by others using 
different methods |2(J . These two distinct scaling regions are obtained according to whether 



the bare persistence length l is larger or smaller than Zosf- The dependence of l\ on l for 
the flexible chains ensures the cross-over, namely the value of the electrostatic persistence 
length for the flexible chain limit crosses over to the Odjik result at l ~ /osf- 

The result for stiff chains (/osf <C lo) is consistent with the assumption that h > k. Thus 
this is valid for wide range of parameters as long as the inequality /osf ^ holds. On the 
other hand, it is not clear if we can neglect the underlined term in Eq. (3.5) for intrinsically 
flexible chains. If we include the contribution from the underlined term, l p may not show 
simple scaling law exhibited in Eq. (3.6). As long as l p ^> l , however, the self-consistent 
equation (3.5) admits a solution for which l p varies as This can be easily seen if we 
rewrite Eq. (3.5) in terms of x = l\ Hi cLS follows 



2 Wo j\ 1 
X ~ — 77T < 1 



A 2 



x 



7T , (\ 

tan -1 - 

2 \x 



l P > lo- (3.7) 



This is a transcendental equation for x which cannot be solved analytically. The solution 
for x, however, is a function of IoIb/A 2 only which does not assume a simple form except for 
a case of l\K 1. In general, we can write 

lp = g {^) /p>>/ ° (3 ' 8) 
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where Q is a function of IqI-q/A 2 only. Thus, the behavior of the persistence length which 
varies as k -1 is the general rule for the case of l p ^> l even though its dependence on l l B /A 2 
is complicated. Explicitly Q satisfies the following equation 



■ 2 Wo f- I 



< 2 gf „ ^ \ 1 



7T , / 1 

tan — - 

2 V2£ 



(3.9) 



which can be solved graphically. As mentioned earlier, the above equation allows a simple 
solution, i.e, Q ~ (IqI-q/ A 2 ) 1 / 2 , if ^> 1 or l p ^> Otherwise we expect a different 

functional behavior of Q. 

According to the above analysis, the persistence length of polyions in the non-asymptotic 
limit, which is assumed to be well described by a trial Gaussian weight with an effective step 
length l\, shows a simple scaling behavior only when l p ^> Iq. On the other hand it is not 
clear if this is valid for other range of parameters. In order to test this we solved the self- 
consistent equation (3.38) numerically for a few cases, with alls/ A 2 ranging from 2-6. The 
values of A, and I were chosen to be unity. All lengths are estimated in unit of I = 1. In 
Fig. 1, In l e is plotted against In k _1 . This figure shows that when al^lj A 2 exceeds 3 then l e 
has a sub-linear dependence on i.e., l e ~ k~ v with y varying from 0.6-0.9. These results 



appear to be in accord with the recent simulations performed by Micka and Kremer [I3|. 
The experimental data of S. Forster et al. also suggest that the persistence length does 
not show a simple power law behavior in certain regimes. When alsl/A 2 is less than about 
3 there is no unique scaling behavior, i.e., l e does not vary as k~ v . 



4 Crossover Conditions and Comparison with Simula- 
tions 

Based on Monte Carlo simulations for weakly charged polyelectrolytes, Micka and (MK) fl3]j 
have arrived at few interesting conclusions. They varied both k and the intrinsic stiffness of 
the chain. Their principle findings are: (a) In the limit of large stiffness, the electrostatic 
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persistence length l e ~ k~ 2 as predicted by the OSF theory. This is valid when the intrinsic 
persistence length exceeds the electrostatic persistence length, i.e., l ^> l e . For flexible 
chains MK find that l e ~ K~ y {y < 1), and y — > 1 as Iq ~ Z e . (2) MK also showed that Z p 
is a strong function of the intrinsic persistence length, l Q . This implies that Iq is altered 
by balance between the bending rigidity and the electrostatic interactions. They suggest 
that the contour distance, s c , at which the contributions from the bending rigidity and 
electrostatic interactions become comparable goes as s c ~ k x with x > 0. 

Our crossover condition for the change in the persistence length goes from the OSF result 
(y — 2) to y — 1 is predicted to occur when Iosf ~ (See Eq. (3.b)). We will argue that 
this is not inconsistent with the findings of MK. According to Eq. (3.6), when Iosf ~ ^0; the 



persistence length l p w l e « V^osf^o ~ This is the crossover condition for the transition 
from y = 1 to y = 2 reported by MK. 

In their paper, MK suggested that the ^-dependence of the contour crossover length 
from our earlier theory gives s c ~ k~ 2 . We believe this is an erroneous interpretation of 



our earlier theory (See Eq. (15) of Ref. [TO]). The crossover length for deviation from OSF 
behavior for l e need not be the same as s c , which describes the contour length scale when 
electrostatic effects become comparable to rigidity. BJ computed s c by considering (6 2 (s)) 
where cos6'(s) = u(s)-u(0). When (9 2 (s)) varies as 1/(Io + Iosf), we expect that the bending 
term is dominant. By identifying the crossover length scale s c to be a length scale where the 
value of (0 2 (s)) crosses over from l/l to l/(l + Iosf), BJ showed that s c ~ 

We can compute s c directly using an variational equation for flexible chains when IqIb/A 2 w 
(9(1), i.e., in the so-called non-asymptotic limit. The contour crossover length can be com- 
puted by analyzing Eq. (3.4). The condition for computing s c is that l e ~ Iq, which implies 
using Eq. (3.4) 

l 18 A 2 Jo J (2vr)3g2 + K 2 e • ^-t) 

This equation is too complicated to solve analytically. However, it is easy to solve Eq. (4.1) 
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numerically. The resulting solution shows that s c is an increasing function of k. The increase 
is not very dramatic which suggests that s c ~ k x with x being small. The parameters used 
in solving Eq. (4.1) are 1 = A = Ib — 1- At least in this case our results are consistent 
with the numerical results of MK. A more detailed comparison between Eq. (4.1) and 
the simulations would require a thorough mapping of the variables in the theory and the 
simulation parameters. We believe that, at least qualitatively, the theory is in accord with 
the computer simulation results. 



5 Conclusions 

The variation of the electrostatic persistence length of flexible charged polyelectrolytes with k 
is complicated. Several regimes can emerge depending on the controlling parameter IqIb/A 2 . 
Here we have inentified the limiting behavior in certain regions of the parameter space. The 
computation of the electrostatic persistence length for flexible chains in the asymptotic limit 
(1 Ib/A 2 <C 1) suggests that the original OSF result is valid. These results were obtained 
by a trial Hamiltonian mimicking a semiflexible chain under tension [TTj together with the 



uniform expansion method [19]. It appears that by focusing directly on the calculation of the 
mean square end-to-end distance effects that are relevant on the length scale on the order of 
persistence length can be implicitly taken into account using Edwards-Singh method. The 



results presented here together with our earlier theory |10[ give a rather complicated pic- 
ture of the dependence of the electrostatic persistence length on the inverse Debye screening 
length. We predict that in the regime when IqIb/A 2 is very large or very small the essential 
features of the OSF result is valid. In the limit of small IqIb/A 2 the blob picture suggested 
by Khokhlov and Khachaturian |22| appears to be valid. In this case, the mean distance 



between the charges is roughly a ^-independent blob size with each blob carrying a renor- 
malized charge. In the non-asymptotic regime the scaling of l e with k is more complicated. 
Here we can expect fractional power for the dependence of l e on k or in some instances 
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complete absence of simple scaling laws. Since Iq and A can be experimentally varied over a 
large regime it may not be possible to analyze the data on the dependence of the persistence 
length of polyelectrolytes using simple scaling laws (Cf. Beer et al. fT8|l ). 
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Appendix A 



In this appendix, we describe the evaluation of the quantities appearing in Eq. (3.18). 
We have argued in the text of the paper that the contributions on large length scale (greater 
than rj) induce an interaction of the excluded volume type, and does not affect the local 
ordering which is described in terms of a persistence length. The average of (exp[z(k + 
q)/r/3 • J* t(s")ds"])' on the short length scale can be approximated as 



In 1 







d 3 q 



A 2 L' J J\s'-s\<^ «9k 2 Jk=o7 ( q -k) 2 + K 2 

« h + h + h 



-(q-k) 2 «| S '- S |/6/ iqir/3- f* t( S ")d S "\' 



where 



In l r s 



h = A-Hrk dS 







Ok 2 



d 3 g 



q 2 + k 



-q 2 ls/6 



I 2 = ^L [ S dss 2 

A 2 r Jo J q 2 + k 2 



q 2 d 3 q c _ q 2 ls/6 



(Al) 

(A.2) 
(A3) 



In 1 



A 2 U J JS<\s'-s\<ri 







Ok 2 



d 3 g 



q 2 + k- 



The first term in Eq. (A.l), I\ can be expanded in powers of k 2 SI (which is much smaller 
than 1) 

d 



1 A 2 I 2 



Ok 2 



i{l-[l-$(K^76)]e^} 



(SA) 2 l B [l + 0{k 2 IS)] 



(A5) 
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where <&(x) is the error function defined by 



2 f X 2 

$(x) = —= e~* dt. 
V 71 " "'o 



(A6) 



Similarly , I 2 , (Eq. (A. 2)) can also be computed. In this case we set k to zero right at 
the outset without any loss of generality. Since this integration involves length scale which 
is much smaller than the interaction between charges within this length scale can be 
assumed to be unscreened. Thus I 2 is approximately given by 

r IB I f S 



\ 2 rL dss2 J d3 ^ q2ls/6 ^(S/A)% 
Finally I 3 can be approximated as 



(A7) 



h = 



A 2 I 2 



dn 2 



— f f dxdx' f 

U J JD<\x'-x\<I p J 



q z + k 1 



rod 





D 2 


- d - 


f 


A 2 


I 2 


dn 2 - 


Jd 



h e 
da;— 



x 



-A' K 2 



(A.S) 



1 + 0{kD) + C(e~ K/ p) 

By following exactly the same procedure the first term in Eq. (3.19) can be approximated 
as sum of two terms, namely, 

d 3 q 



Ib I 2 



( -^-q 2 (s> - S ) V^l^l/ 6 (e^ 3 -r t ^ ds ")' = h + h (A9) 



A 2 L' J JS<\s'-s\<ri J q 2 + K 2 

where J 4 is the contribution to Eq. (A. 9) coming from the scale \s' — s\ < S and is given by 

S 2 



A 2 ' 



(A10) 



The term J 5 arises from contributions on the large length scale S < \s' — s\ < rj and can be 
approximated as followed 



Ib D a 1 



A 2 I 2 U J Jd<\x'-x\<i p 
A 2 I 2 



q 2 + k 2 



rod 



l B rh 2 d 2 e 
/ axx 

JD 



dx 2 x 
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l B S 2 /A 2 (All) 



APENDIX B 

To illustrate the limitation of the Gibbs-Bogoliubov variational approach in the compu- 
tation of persistence length, let us consider a simplified system which is described by the 
following Hamiltonian 



TL Iq 



£ 9 2 (s)ds + \£ £ 0(s)K(s, s')9(s')dsds' (B.l) 



k R T 2 



where cos9(s) = u(s) • u(0) and 9(s) = The kernel K(s,s') describes interaction 

between charges on the chain and has yet to be specified. If the kernel K(s, s') is 

S2S! (S2-Sl) 



A 2 Jo Js+k ds 



then the Hamiltonian in Eq. (B.l) corresponds to a polyelectrolyte chain near the rod limit. 
For the present purpose, we use a simpler form of K(s, s') whose Fourier transform is given 



by K q = ^flosF- 



H 



k B T 2 „ 

Since the Hamiltonian in Eq. (B.3) is quadratic in the variables 9 q all quantities with refer- 
ence to Eq. (B.3) may be exactly computed. The efficacy of the approximate methods can 
be assessed by comparing to the exact result. In the following, we adopt three different trail 
Hamiltonians in the computation of the persistence length to illustrate possible limitations 
in each case. 

(a) First, let us use a trial Hamiltonian 

Tit 



= ±E? 2 >W- S . (5.4) 
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Within the Gibbs-Bogoliubov variational approach, the best choice for K q is obtained by 
minimizing the variational free energy T v = T t + (H — Ti-t) t - Performing the Gaussian 
integration over 9(s), we have 

Minimization of the variational free energy with respect to K q leads to 

*. = «. + fe»VW (B - 6) 

When a general form of variational kernel is used, the variational model produces the same 
result as the original Hamiltonian. 

(b) As a second example, let us choose a bending model as a trial system which is 
described by 

z q 

This is such that (0 2 (s)) = s/l p . In this case, the variational free energy is given by, neglecting 
numerical prefactors, 

^ ~tt + T"Iosf (B.8) 



k B T Al p l p 

where A -1 is the upper bound for q. Minimization of the variational free energy results in 

ip ~ io + (B.9) 

When a bending model is used, the Gibbs-Bogoliubov variational approach leads to a wrong 
result for the electrostatic persistence length of a polyelectrolyte chain near the rod-limit 
which scales as kT 1 . This is because the free energy of the bending model cannot be made 
mimic the free energy on scales on the order of the persistence length. In other words, 
short length scale (l p ) properties in the original system are washed out by minimizing the 
variational free energy with respect to l p . 

(c) Let us now compute the persistence length of a polyelectrolyte chain described by H 
in Eq. (B.l) using the Edward- Singh's self-consistent method. The following derivation is 
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analogous to that given in the text (Cf. Eq. (3.11)). To this end, we write 

and compute (9 2 (L)) self-consistently by considering the underlined term to be small. Up 
to the first order in the underlined term, the bending model can be made the same as the 
real system by requiring 

(^o-g|-^ 2 W) t = (^(^)|E^X^) t -^ 2 W) t (lE^X^)r (B.ll) 

(Up q q 

A little algebra leads to the following equality from the Eq. (B.ll) 

uip i p q q 

For large L most of contribution to the sum on r.h.s of Eq. (B.12) is from q ~ 0. This results 
in 

l P ~l + losF- (5.13) 

Thus the ES method provides a reliable way of calculating the persistence length l p . The 
result using this method agrees with the exact answer for the simple Hamiltonians of the 
form given in Eq. (B.3). By focusing directly on the local quantity the ES method implicitly 
includes effects on the scale of the persistence length. Such fluctuations are integrated 
out when the free energy is computed. Consequently, the persistence length needs not be 
accurately given by minimization of an approximate free energy. 
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Figure Caption 



Fig. (1). Plot of the electrostatic persistence length, / e , as a function of k~ 1 . Both quantities 
are expressed in natural logarithm values. The various curves correspond to different values 
of the parameter alsl/A 2 (See Eq. (3.5)). The values of a from the top to the bottom curves 
are 6, 4, 3, 2.5, and 2 respectively. For alsl/A 2 > 2.5, l e oc k~ v and the exponent can be 
obtained from the slope of the displayed curves. The exponent y varies from 0.6-0.9. For 
cxIbI/A 2 < 2.5, there appears to be no simple power law describing the dependence of l e on 

K. 
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